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Abstract 

It is postulated in general relativity that the matter energy-momentum tensor 
vanishes if and only if all the matter fields vanish. In classical Lagrangian field 

7-H . theory the energy and momentum density are described by the variational (sym- 

metric) energy-momentum tensor (named the stress tensor) and a priori it might 
occur that for some systems the tensor is identically zero for all field configurations 
whereas evolution of the system is subject to deterministic Lagrange equations of 

*q | motion. Such a system would not generate its own gravitational field. To check 

if these systems can exist in the framework of classical field theory we find a rela- 
tionship between the stress tensor and the Euler operator (i.e. the Lagrange field 



equations). We prove that if a system of interacting scalar fields (the number of 
fields cannot exceed the spacetime dimension d) or a single vector field (in space- 
times with d even) has the stress tensor such that its divergence is identically zero 
(i.e. "on and off shell"), then the Lagrange equations of motion hold identically 
too. These systems have then no propagation equations at all and should be re- 
garded as unphysical. Thus nontrivial field equations require the stress tensor be 
nontrivial too. This relationship between vanishing (of divergence) of the stress 
tensor and of the Euler operator breaks down if the number of fields is greater than 
d. We show on concrete examples that a system of n > d interacting scalars or 
two interacting vector fields can have the stress tensor equal identically zero while 
their propagation equations are nontrivial. This means that non-self-gravitating 
(and yet detectable) field systems are in principle admissible. Their equations of 
motion are, however, in some sense degenerate. 

We also show, that for a system of arbitrary number of interacting scalar fields 
or for a single vector field (in some specific spacetimes in the latter case), if the 
stress tensor is not identically zero, then it cannot vanish for all solutions. There 
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do exist solutions with nonzero energy density and the system back-reacts on the 
spacetime. 

PACS numbers: 04.20.Fy, ll.10.Ef 

1 Introduction 

It is commonly accepted in general relativity that the gravitational field is generated 
by all forms of matter (i.e. all species of elementary particles and fields) and that all 
features of a given form of matter that are relevant for determining its gravitational 
field are encoded in its energy-momentum tensor (cf. e.g. |H 12])- In other terms en- 
ergy and linear momentum are the source of the gravitational field. This postulate 
does not tell one how to construct the energy-momentum tensor for a given kind of 
matter and whether it is unique. In principle this tensor should be determined in a 
special-relativistic theory describing the material system and then minimally coupled 
to gravity. However a reliable expression for this tensor has been found only in few 
cases (classical electrodynamics, relativistic hydrodynamics etc.) and it has turned out 
that the canonical energy-momentum tensor, though being conceptually important due 
to the first Noether theorem, does not provide in most cases the correct value of en- 
ergy density and therefore is unphysical. A unique universal definition of energy and 
momentum density arises if the equations of motion for the matter under consideration 
can be derived from a Lagrangian. This is the Hilbert's variational (with respect to the 
spacetime metric) energy-momentum tensor (or Belinfante tensor), hereafter denoted 
as the stress tensor. It is worth emphasizing that the physical energy-momentum tensor 
cannot be determined by merely manipulating with the matter equations of motion and 
the use of the stress tensor is indispensable [3] . Now it is commonly accepted that it is 
this tensor that correctly describes energy and momentum density and their flows for 
any Lagrangian matter both in curved and flat spacetimes jl]. 

It is postulated that the Lagrangian for any matter is such that the resulting stress 
tensor T^ v vanishes on an open domain in the spacetime if and only if the matter fields 
vanish on the domain |T] . This condition expresses the principle that any matter carries 
energy. The energy conditions ([1J, Chap. 4) actually imposed on a matter Lagrangian 
exclude negative energies and ensure the "only if condition. 

In this paper we investigate the problem whether there is some classical matter hav- 
ing no energy at all, i.e. whether there exists a matter Lagrangian giving rise to the 
stress tensor vanishing identically, T^ u = 0, for all field configurations independently of 
the field equations ("on and off shell"). Such a matter should be subject to deterministic 
equations of motion (the Lagrange ones) and according to Einstein field equations would 
propagate as a test one in a fixed spacetime. This matter would be non-self-gravitating 
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and its interactions with other fields would be severely restricted by the constraint that 
they should exclude any energy transfer. At first sight one might conclude that the 
lack of any energy exchange with ordinary matter systems implies that this matter is 
nondetectable and as such it may be merely ignored. This is the case in quantum theory. 
However, here we are dealing with classical fields and in classical physics any system 
may be treated as open, i.e. as being in contact with some surrounding. The contact 
may be arbitrarily weak, i.e. involve negligibly small energy transfer (or no transfer at 
all), nevertheless by observing the surrounding (which acts as a "marker") one can make 
measurements on the system. In Example 1 we give a hint of how a system of scalar 
fields without energy might be detected by their coupling to the electromagnetic field. 
We do not pursue the problem further, it is sufficient to say here that the lack of energy 
does not imply "physical nonexistence" in the sense of nondetectability. 
This is why we investigate in this work to what extent Lagrangian field theory admits 
such bizarre systems and whether some of them can be excluded on theoretical grounds. 
To this end we show for some classical fields that there is a close connection between 
energy (the stress tensor) and deterministic equations of motion. For these fields we 
prove that the identically vanishing stress tensor implies that the Lagrange equations 
also identically vanish for all field configurations, i.e. the fields are subject to no prop- 
agation equations at all. On very generic physical grounds one can then reject these 
fields as unphysical. 

We first formulate the problem in full generality. Consider a system of classical 
matter tensor fields if) a with a collective index A (there is no need to deal with spinor 
fields since they can be expressed as tensor ones An example of a vanishing stress 
tensor is provided by a Lagrangian for if) a of the form 

L&A, i>A;n) = -^=L (iff A , ifj A ,»), (1) 



where Lq is a scalar density of the weight +1 (so that L is a genuine scalar) and L$ is 
independent of the spacetime metric and its first derivatives, 

OL dL =Q (2) 



(Here ipA;n = V^Va is the covariant derivative w.r.t. g^ u ). Then the stress tensor (the 
signature is — h ++) 

TMa) = ^=^(V=£) 0) 

vanishes identically and the spacetime evolves, in absence of other forms of matter, as 
an empty one, = 0. We shall consider in the next section a few specific forms of Lq 
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for scalar and vector fields and find in which cases the corresponding Lagrange equa- 
tions are trivial, making the fields unphysical. However it is clear that the conditions 
(flj) and (2) are merely a (restrictive) sufficient condition for having T^ v = and not 
a necessary one. In fact, there are altogether 50 conditions (2) whereas there are only 
10 identities T^i/ja) = actually imposed on a possible Lagrangian. In what follows, 
we shall assume in Propositions 1 to 4 that the field Lagrangian L^a^a-,^, 9^u) m ay 
depend on the metric both explicitly and implicitly (via the covariant derivatives) and 
only in Examples 1 to 5 we shall consider various Lagrangians satisfying eqs. (1) and (2). 

To avoid any confusion we emphasize that we work in the framework of classical field 
theory and thus we do not take into account the classical topological field theories 
such as BF theory .6J. They are metric independent and in this sense they might seem 
relevant to the present work, but their actions are typically given by surface integrals 
and they describe global observables related to the topological invariants of the man- 
ifolds on which they are defined. These theories have no local degrees of freedom, so 
there are no propagating field excitations (particles). 

One may expect that if the system under consideration consists of a large number of 
fields, then vanishing of the full stress tensor will turn out to be a condition too weak to 
make trivial the large number of (coupled) Lagrange equations for the fields. We shall 
show on concrete examples that this is the case. Thus, we should theoretically allow for 
some specific systems of fields which carry no energy nor momentum, which nonetheless 
obey deterministic equations of motion (causal or not). The purpose of the present work 
is to show how many fields and of what type are necessary to this aim and how peculiar 
their Lagrangians must be. 

The main thrust of the paper are Propositions 1 and 2 in sect. 2 and Propositions 3 
and 4 of sect. 3. Propositions 1 and 2 state that if a classical system consists either of a 
number of scalar fields or a single vector field and the stress tensor for the system is zero 
on and off shell, then Lagrange equations hold identically. This system is thus regarded 
as unphysical. However, if the number of scalar fields exceeds the dimensionality of the 
spacetime or there are two (or more) vector fields, the theorems break down. It is also 
relevant whether the dimensionality is even or odd. We show on specific examples that 
for sufficiently large number of fields (or their components) there exist Lagrangians sat- 
isfying conditions and (J2J) and possessing the Lagrange equations of motion. Their 
dynamics is however quite bizarre: for systems of interacting fields there are no free-field 
solutions, for some cases the systems of propagation equations are either degenerate (of 
first order) or indeterministic (less equations than degrees of freedom). One concludes 
that classical fields without energy, though not excluded by principles of Lagrange field 
theory, require very peculiar and easily recognizable Lagrangians and these are unlikely 
to appear in modelling the physical reality. 
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Propositions 3 and 4 of sect. 3 solve, for interacting scalar fields or a single vector 
field, the problem of whether the stress tensor may vanish for all solutions while it is 
nonzero for some field configurations off shell. It turns out that there are always some 
solutions for which the stress tensor cannot vanish. In this sense nontrivial equations of 
motion imply in most cases a back reaction of the system on the spacetime. Unfortu- 
nately, for a vector field the proof works only in a small neighbourhood of flat spacetime 
and can be generalized solely to spacetimes admitting a covariantly constant Killing 
vector. 

In section 4 we make some critical comments on the canonical energy-momentum ten- 
sor. Proofs of Propositions 1 and 2 employ the second Noether theorem and Proposition 
4 requires the Belinfante-Rosenfeld identity. A detailed derivation of both the identities 
is provided in the Appendix. 

2 Scalar and vector fields having no energy 

According to Introduction one should not expect to eliminate as unphysical the systems 
carrying no energy if they involve many interacting scalar and vector fields or fields with 
high spins. Instead one should separately investigate systems involving rather a small 
number of scalar fields or a single vector field. These fields are defined on a curved 
spacetime and are viewed either as test fields or as a part of a larger matter source of 
gravity. No specific gravitational field equations are assumed. The two systems are dealt 
with in the following two propositions. Let d > 3 be the dimensionality of the spacetime. 

Proposition 1 

Let a material system consist of n < d interacting scalar fields and let the divergence 
of the stress tensor vanish identically for all values of the fields, T^ v '. v = 0. Then their 
Lagrange equations of motion also hold identically for all values of the fields. 

Proposition 2 

Let the dimension d be even and let a material system consist of a single vector field 
with a Lagrangian which can be expanded in a Taylor series in At^j in the function 
space of all antisymmetric tensor fields F^ u = —2A\p. v ] defined on an open domain in 
the spacetime. (The series is centered at F^ u = and thus the Lagrangian and all its 
derivatives are regular at this point.) If the divergence of the stress tensor vanishes 
identically for all values of the vector field, T^ v '.„ = 0, then the equations of motion hold 
identically. Moreover, in d = 4 the stress tensor itself is zero, T IXV {A) = 0. 

Idea of the proof. 

The proof of both Propositions is based on the Noether identity, valid for any matter 
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field, which arises from the coordinate invariance of the matter action integral. The 
identity is derived in the Appendix. It involves the divergence T^ v rather than the 
stress tensor itself and this is why the assumption of Propositions is apparently weaker 
than Tjjj, = 0. At least in the case of a vector field and d = 4 the two assumptions are 
equivalent. 

I. Proof of Proposition 1 

Consider a system of n interacting scalar fields <f) a , a = l,...,n. For scalars the co- 
efficients Za^o introduced in Appendix, eq. (A. 15), are zero and the Noether identity 
(A. 17) reduces to 

E a <\> a ^ = (4) 

where depends on all the scalars and their first and second (in the case of a nonmin- 
imal coupling) order derivatives and the same holds for the n scalar quantities E a . (If 
T^v does involve a; /^, then the terms containing third order derivatives, arising on the 
r.h.s. of (4), do cancel each other.) For n < d and arbitrary fields (f) a , the vectors a;/J 
are linearly independent. On the other hand the condition T Atl / ly = gives E a $ a .^ = 
suggesting that the gradients are actually linearly dependent. The consistency is re- 
stored only if E a = 0, i.e. the field equations hold trivially. 

The theorem breaks down for n > d as is seen from the following example 1 . (Our 
experience with some readers shows that it should be explicitly stated that the deter- 
minant of any mixed tensor Y^ v is an absolute scalar, hence |det(Y^)| 1//2 is a scalar 
density, i.e. transforms as y/—g. Accordingly, the Lagrangians in the Examples 1, 4 
and 5 are absolute scalars. On the other hand the Lagrangians in the Examples 2 and 
3 are chosen as pseudoscalars (i.e. they transform as scalars multiplied by J/\J\, where 
J is the Jacobian of a coordinate transformation) merely for computational simplicity 
and can be made scalars by taking the absolute value; this change will not affect the 
conclusions which follow from them.) 

Example 1. 

For a system of n scalar fields one defines = J2a=i 4>a,iJ, ( l > a,v For arbitrary scalars 
and n > d its determinant det(P AtJy ) ^ and may be used to make up a Lagrangian, 
specifically, 

L (/) = - 7 L|det(P Ml/ )| 1/2 . (5) 

Then the conditions ©-(j2J hold and T^ v = 0. Let d = 4 for simplicity, 
i) For n = 4 the Lagrangian takes on a simpler form, 

L 4> = 7y= |det(0 OiM )| 
1 Examples 1 and 4 were suggested to the author by Andrzej Staruszkiewicz. 
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and furthermore it is a full divergence, 



L<f> = V Q (e Q ^>i0 2 , /3 3 , At 04,,) , 



(6) 



where the antisymmetric Levi-Civita pseudotensor (i.e. it transforms as a tensor 

times J/\J\) with £ 0123 = Clearly then E a (<j> h ) = 0. 



ii) Yet in the case of n = 5 scalars the Lagrangian (5) cannot be simplified to an 
analogous form and furthermore it is not a divergence. The five Euler operators do not 
vanish, 



where Q a @ is the symmetric inverse of P^, Q m P av = $„■ To show that E a are not 
identically zero one finds counterexamples. In Minkowski spacetime one puts four of <j) a 
equal to the Cartesian coordinates and the fifth scalar equal to some nonlinear functions 
of time, e.g. t 2 or e*. Then two of the operators are different from zero. One concludes 
that the system of 5 scalar fields with the Lagrangian © is a form of non-self-gravitating 
and carrying no energy matter, subject to nonlinear propagation equations. The rela- 
tionship E a cf) a ^ = shows that there are at most four independent equations for the 
five scalars, thus the equations of motion and the initial data do not uniquely determine 
the evolution of the system. 

The system might be detected, at least in principle, by its influence on the electromag- 
netic field A^. Let j M = qJ2a=iQ fJ ' U( l ) a,u, where q is a coupling constant, be a current 
associated with the scalars. The current is coupled to the electromagnetic potential via 
the standard interaction term, Lj n ^ = j^A^, and the stress tensor corresponding to the 
Lagrangian + Lj n ^ remains identically zero. Yet the electromagnetic field is affected 
since the current enters the Maxwell equations. 

II. Proof of Proposition 2 

Here the main idea of the proof consists in expanding the Noether identity for the vector 
field, in the case where T^ v ;u = 0, in a series of identities (not equations) ultimately 
resulting in the Euler operator which vanishes identically. 

For a single vector field the Noether identity (A. 17) takes the form (A. 18), 
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(7) 



6=1 



A a - t p being a "field strength" and the Euler operator 



(8) 



with F a p = Ap.a - 
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We set T au ' u = and decompose the identity (8) into a sum of terms containing covariant 
derivatives of of definite order. To this aim we first introduce tensors 

r) 2 T 

f>H~/af3 — w _ p>a/3fj,-y (->r>\ 

d ~dA a ,, and 1 -dAS (U) 
these are functions of A^ and A^. v (and the metic). Then one finds 

Bd^ a 

E , = - B ^A a , M - -^A v , a + F. (12) 

We insert this expression for into (8) and decrease the order of derivatives of A^ with 
the aid of Ricci identity. After some manipulations the identity E^F ail — A a E^.^ = 
takes on the following involved form, 

= A a B^ u 1 Ap. l{ ^ v) + A a -gj—A x . {Ma) A u . m + 



A\;(tur) { -F av B uaX ^ + A 



A 7 Rru/3-y — — h 2 — h 



dA V]P dA p dA\ dA, 

-\-a a (A^, Ap-v, Rf iu \a)i (13) 

where complicated term containing no higher derivatives than of first order. 

Being an identity for arbitrary A^, the terms with derivatives of different order cannot 
cancel each other, instead they should vanish separately. This implies that the identity 
splits into a cascade of four sets of independent identities. We first consider the term 
with third derivatives. It vanishes if their coefficients are zero. Since A a ^ one gets 
B^fh = 0, then the definition (10) yields 

These in turn imply a specific relationship between B a P^ u and d a ^. It turns out useful to 
decompose A a .p and d a/3 into symmetric and antisymmetric parts, A a .p = A( a; ^ + \F$ a 
and 

( fi°n = 0£— = S°t l , cM=2^- = N^. (15) 
oA^.p) oFp a 

Since 



8A^ U dA a;/3 
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one finds 

dS«P dN<# dST* dN*# 

dA M ^dA M dF^ dF, u 1 ] 

and a similar expression with the pairs af3 and jj,u interchanged. By antisymmetrizing 
in these pairs one gets 

sMH = _ 2 dN^ = (17) 

On the other hand the symmetries (14) mean that B a/3fMU = .B^HH and by equating 
eq. (16) to (17) one arrives at the following restrictions imposed on <i Q/3 , 

Qgap Q S aP Q N aP ^ 



From these one infers that the Lagrangian may depend on A^ a .^ only via linear terms, 

L = his^A^A^ + Hs)g al3 A (a;P) + L 2 (A„ F^), (19) 

where f\ and fa are arbitrary smooth functions of the vector length, s = A^A^. The 
first two terms in (19) can be expressed as V a (f(s)A a ) + h(s)A? a , with 2d//ds = f\ 
and h = ji — f. Then discarding the divergence term from the Lagrangian one finds 

8N a/3 8N^ V 

d a ? = hg ap + N ap and B a ^ v = -2—— = -2——. (20) 

oF^ oF a p 

We next study the identities involving A v .^^ quadratically. From (13) one sees that 
vanishing of the coefficients of these terms requires 

^T + ^T s0 ' (21) 

By employing eqs. (20), however, one cannot simplify these identities (after inserting 
(20) into (21) one recovers the identities with the pairs of indices ^7 and pia inter- 
changed), therefore one passes to studying the identities involving Ax-^a) linearly. These 
imply vanishing of the curly bracket in eq. (13). Making use of (20) this reads 

F av B v ^ x + A a C^ x = (22) 

with 

C^ x = 2ti(g^A x - A^g a)X ) - - ^— + ^— = C^ a)x (23) 

2 y oA^ oA a J 



and h! = dh/ds. All the time one investigates a generic field for which det(-F)^) 7^ 0. 
(It is here that the even number of spacetime dimensions becomes relevant; for d odd, 
det(F flu ) = 0.) Then there exists the inverse matrix f a/3 = —f^ given by f^F^p = 
F^f = 5%. Multiplying (22) by f ar one gets 

One immediately sees that this expression implies A a B a ^ u ^ = 0. On the other hand 
one infers from (20) that B a ^ u ^ = B^ 3 ^^ . Applying this symmetry to the former 
identity one finds 

B <*Qw)P A/} = -f^AvC^Ap = 0. (25) 

This in turn implies that C^Ap = since for a generic vector field f aa A a 7^ 0. By 
inspection of the expression (23) one concludes that the term proportional to h! in the 
latter identity contains no derivatives, whereas the other term in this identity is a sum 
of terms each of which does contain F^ u . In fact, terms of the form dN^ v / dA a cannot 
involve an additive term free of the derivatives, since the latter would only arise from 
a linear term in the Lagrangian, k^ u (A)F^ u . However, for d even, an antisymmetric 
W cannot be made up alone of A^, and the Levi-Civita tensor. In conclusion, the 
identity C^^Ap = splits into two sets, 

h'ig^s - A^A u) ) = (26) 



md Mw + im7H (27) 

Eq. (26) may be satisfied only if h — ho — const; then hoA? a is discarded as being a full 
divergence and finally L = L(A^, F^ u , g^ u ). 

We now return to investigating the expression (24). It is very peculiar. First, the 
r.h.s. of the expression does not have the above mentioned symmetry B a ^ u ^ = B^ v ^ a . 
Second, according to the assumption of the theorem, B atJiUf3 is analytic (is a Taylor series) 
in F^ about F^ = and the same holds for . Yet B a{ -^ is proportional to f aa A a 
and the expression should be valid also about F^ u = where the inverse f aa does not 
exist. (Actually the expression breaks down whenever det(-F M; ,) = 0.) The dependence 
on f aa cannot be eliminated since there is no contraction of this tensor with any of F^ 
appearing in C^ v<i . One then infers that identity (24) may hold only if B a ^ u ^ = 0. 
To avoid any confusion it should be emphasized that the above reasoning does not apply 
to the Noether identity (8). In fact, with the aid of f afi it can be reexpressed as 

E a = n(A^ v +T, u n 
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and apparently E a does explicitly depend on f afl , contrary to the assumption of the 
theorem. However, here f a>1 is contracted with T^ v and the original form of the iden- 
tity ensures that the /"^-dependence is trivially cancelled. This fact stresses the role 
played by T^. 

The identity B a ^ u ^ = makes identities (21) trivial and has two further conse- 
quences. First, together with (14) and (20) it implies that B a ^ uf3 is totally antisymmet- 
ric, 

£>anu/3 _ g[anvl3\ ^g) 

In d — 4 one then infers that B afil/ P = p(A\, F\ a )e afiUl3 with some definite pseudoscalar 
function p. Second, it requires = or from (23), 

dN a » dN av . . 

aI7 + W"°' (29) 

The latter shows that the following tensor is totally antisymmetric, 

dN afl r i 

n^E — =nK (30) 
a A v 

The formula (12) for the Euler operator is reduced, upon applying eqs. (28) and (30), 
to 

E» = -^n^F a p + F, (31) 

now it contains no second order derivatives. 

Finally we study the last system of identities arising from (13), those involving at 
most the first derivatives of A^, i.e. a a = 0. After some manipulations with the use 
of the Riemann tensor symmetries and eqs. (28) and (30) one arrives at the following 
expression: 

dE* 

a a = F aik E» - A a g^-A^ ee (32) 

with E^ given by (31). As in the case of one may use the analyticity property to 

prove that E^ = 0. It is interesting, however, to see that this result can also be attained 
in an independent way. Using the formula (31) one easily finds that the scalar in the 
last term of (32) is equal to 

dE^ d 2 L A 

dA u Av ^~ dA,dA u A ^- {dd) 
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One sees that the r.h.s. of eq. (33) depends linearly on A^.^ and upon inserting the 
scalar back into eq. (32) the term F ail E p acquires the same dependence. On the other 
hand it has already been proved that L = L(A^, F pu , g pu ) and the symmetrized derivative 
cannot arise in the process of differentiation of the Lagrangian. The contradiction is 
removed by requiring that 

a%k s0 (34) 

or L = L^{F a ^)A p + L (F a p). However for even number of dimensions it is impossible 
to make up a vector out of F a p, g pu and the Levi-Civita tensor. In consequence L M = 
and the Lagrangian is some function of F a p and alone. 

As a result of vanishing of the scalar (33) one gets from (32) that F a/J ,E^ = and finally 
one arrives at the conclusion that for a generic vector field the Euler operator associated 
with the Lagrangian generating the stress tensor with T pu ]V = vanishes identically, 
E»[L(A)] = 0. 

This outcome is in agreement with eq. (31) since the Lagrangian does not depend on 
A^ and then n * 4 " = = Z*». 

The last step of the proof consists in finding a relationship between the identity 
T^ u ' u = and the stress tensor itself in the case d — 4. From eq. (28) one infers that 

o2 T 

B a ^ = 4 — — = p(F Xa )e a ^ (35) 

with unknown pseudoscalar p. The Lagrangian, which depends on A a only via F a p, is 
a function of the two invariants of the field strength, V = F a/3 F a ^ and W = P 2 where 
P = e^p^F^ is a pseudoscalar, i.e. L = L(W, V). Then 

B a»v/3 = 32PL wv (e a < lX(7 F u > 3 + F a "e^ Xa )F Xa + 32{L W + 2WL ww )e apXa F Xa e ul3Tp F rp 
+16L vv F a ^F v(3 + AL v (g au g^ - g af5 g pu ) + \QPL w e a ^ 
= J a ^ + 16PL w e a ^ = p{W, V)e a ^ ul3 , (36) 

where Ly = dL/dV etc. The first four terms have lower symmetry than that required: 
only ajj,uf3 = [a/j] [vf3] = i/fla/i, thus their sum J a vP must vanish identically and this is 
possible only if each of them is separately zero. To see this one assumes that the sum is 
zero for a fixed field F 0pu . Let F pv = F 0pu + 5F pu where 5F pv is arbitrary infinitesimal. 
The variation 5J aflul3 also must vanish, 

sja ^p = 325(PL wv )(e a » Xa F^ + F^e^ Xa )F 0Xa + 

32PL wv [{e a ^ X(J 5F^ + 5F ap e ul3X(7 )F 0XcT + {e apXa Ff + F^e upXa )5F Xa ] 
+325(L W + 2WL ww )e a ^F 0Xa e^ Tp F 0rp 

+32(L W + 2WL ww )e a ^e^ Tp (5F Xa F 0Tp + F 0Xa SF Tp ) + 1Q5L VV F^F^ 
+ 16L vv (5F ap F^ + F^5F ufS ) + A5L v (g au g pl3 - g^g^) = 0. (37) 
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Here for any $(W, V) one has as usual = $ W 5W + $ V 5V where 5W = 2P 5P = 
4P E al3 » u F 0al3 5F^ and 5V = 2Fg u 5F^ u , all the derivatives and W and V are taken 
at Fq^. At fixed Fq^ one varies SF^ in such a way that the variations 5W and 8V 
remain unaltered. Then the first, third, fifth and seventh (the last) terms in (37) remain 
constant whereas the other three are variable. For the identity holds for any SF^, the 
constant terms must vanish, i.e. 5(PL WV ) = 5(L W + 2WL WW ) = 5L VV = 8L V = 
and these identities hold for any values of 5W and 5V. The last identity yields then 
Lyw = Lyv = or L(W, V) = L(W) + aV, a is constant. Inserting this Lagrangian 
into eq. (36) one sees that the first and third terms are zero whereas the fourth term 
becomes independent of F^, thus it must vanish too. This implies Ly = or a = 0. 
The vanishing sum J a v v P = q consists now of the second term alone and it is zero 
provided L w + 2WL WW = 0. A general solution to this equation is (dropping an 
additive constant) L = cyW = c\P\ with constant c. Finally, 

L = 2c\V a (e a ^A^)\, 

a textbook result [7j. 

This completes the proof of Proposition 2. 

For dimensions d = 2n, n > 2, the tensor F a p has more invariants. One may 
conjecture that also in higher dimensions eq. (28) implies that 

r „\ c ai/3i...a n /3 n Z7 T? I 

L - C \ £ *aiPi ■ ■ -^a n p n \ 

being a full divergence and hence T liV {A) = 0. However it is harder to prove that this 
solution is unique. 

In principle one might envisage a scalar (or n < d scalars) or a single vector field in 
d > 4 even with ^ and T M1/ :Jy = 0. Formally such a field might appear in Einstein 
field equations. However, since the Lagrange equations of motion are trivial and the field 
has no determined propagation, it should be rejected on physical grounds. At first sight 
the stress tensor for a ground state solution of a quantum field, < T^ w >= pvg^u, where 
Pv is a constant energy density for the classical nonzero value of the quantum field in 
this state, contradicts the above statement. It should be therefore emphasized that this 
expectation value of T^ u comes from semiclassical considerations and cannot be derived 
within classical field theory (formally this expression is generated by the cosmological 
constant term in the full Lagrangian including gravity); the proposition does not apply 
to that case. 

The Proposition 2 cannot be generalized neither to odd number of dimensions nor 
to more than one vector field. This is shown by the following counterexamples. 
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Example 2. 

Let d = 3. If one chooses a pseudoscalar L = e a ^F a pA J: then T^A) = 0. This 
Lagrangian is not a divergence and = 2e^ a/3 F af 3, the operator involves no second 
order derivatives since L is degenerate being linear in F a p. The field equations read 
then F a/ 3 = and admit only one solution, the vacuum. The model is trivial. 



Example 3. 

In d = 5 the analogous model is nontrivial. Let 

L = W^A, = e^ 5 F a/3 F lS A„ (38) 

where = e m ^ lS F^F^- The stress tensor is zero, T^ V {A) = 0, whereas the La- 
grangian is not a divergence and is gauge invariant under A'^ = A^ + d^x- L is a P~ 
parently nondegenerate as being quadratic in F a p, nevertheless the Euler operator is 
of first order, E^ = 3W^, and in this sense the propagation equations for the field are 
degenerate. One finds that Eft = owing to the gauge invariance. Then the Noether 
identity (8) reduces to F^E^ = with E^ nonvanishing in general (det(F^) = 0). The 
field equations = are quadratic in F a p and admit nonzero solutions. A particular 
solution is, e.g. 

A l = axx - a 2 x 2 - a 3 x 3 - a 4 x 4 , A = A 2 = A 3 = A 4 = 0, (39) 

yielding F i = a±, F\ 2 = a 2 , -F13 = a 3 and Fu = a 4 , otherwise zero, with a±,...,a4 
constant. 



Example 4- 

Next we consider two vector fields forming an open system, d = 4. A field A^ in- 
teracts with a given external field W^. First one defines an antisymmetric tensor 
= W a A [a d^A v] = \W a A [a F^ v] , where as usual, F^ = A v ^ - A^ v . The tensor 
is metric independent and thus may be used to constructing an interaction Lagrangian, 

L(A,W) = -L[- det(^)] 1 / 2 = ^-[2V a pV^V» v V™ - (^F^) 2 ] 1/2 
V ~~ 9 V 8 

= \\s a ^ v V aP V, v \. (40) 

L does not involve derivatives of the external field. In this model there is no free 
Lagrangian for A^ and obviously L(A, 0) = 0. The stress tensor generated by the 
interaction Lagrangian is T^ U (A, W) = 0. On the other hand if one attempts to complete 
the system, i.e. to make it closed, by addding a Lagrangian L w for (free or including 
an interaction term), then y/—gL w must depend on the metric. Thus a full stress tensor 
for a closed system of the two vector fields is different from zero. 

The open system described by L as in (40) has no simple symmetries. The Euler operator 
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for is nondegenerate and very complicated. The second order derivatives appear 
in it in a term of the form 



e^A v A a . m WtW°A a . (41) 

Example 5. 

One can construct a model for a closed system of two interacting vector fields in d = 4. 
One assigns field strengths to vector potentials A^ and B^. 

Fjxv — — A^ v and = B v<fl - B^ v . (42) 

The interaction Lagrangian is chosen as 

L(A,B) = -L[det(F /1J/ )det(^)] 1 /4 = \\L A {F)L B {H)\V\ (43) 

where La = e^^F^F^ and L B = e af)txv H^H^ are pseudoscalars. As in Example 4 
there are no free Lagrangians for the fields (by introducing free Lagrangians one will 
render the full Lagrangian density metric dependent, Proposition 2) and the interaction 
Lagrangian is gauge invariant under A'^ = A^ + d^Xi an d B' = B^ + d^X2, Xi an d X2 
arbitrary. The interaction energy and momentum for the system are zero, T^ U (A, B) = 0, 
whereas the Euler operator consists of two vector ones, 




and analogously, by applying the symmetry of interchanging the fields in L, 

ft - -Zr" 1 - (tt - L -t) ■ (45) 

The operators contain second order derivatives (are nondegenerate). One can find a 
special solution to the field equations E A = = Eg in Minkowski spacetime in the 
following way. One assumes that the "electric" and "magnetic" fields corresponding to 
F^ v are parallel and of equal length and direction, then F^ V F^ V = and L A = —8a 2 , 
a > with Fqi = F 32 = a, otherwise zero; and analogously, H^H^ = and Lb = —8b 2 , 
b > and H m = H 32 = b. In this case the two operators are proportional, 

K = \{-d u do, d 3 , -d 2 ) In a - = - a -E%. (46) 

A particular solution is then A^ = —a(x,0,0,y) and B^ = —b(x,0,0,y) with a and b 
constant. Eq. (46) shows that a and b cannot vanish (no free fields). 
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3 The stress tensor vanishing for solutions 



The Propositions 1 and 2 raise an obvious problem: is it possible that T^, while being 
different from zero for some values ("off shell") of the field under consideration, if) a, 
vanishes for all solutions of E A = and not only for the trivial solution ipA = 0? Such a 
possibility would be more disturbing than the case T^ u = for all values of the field, for 
field configurations which do not obey some equations of motion do not exist in nature. 
Such a matter, subject to deterministic causal propagation equations, would be a truly 
nongravitating one. In the framework of classical field theory one can almost exclude 
such fields under reasonable assumptions. We again consider only scalar fields and a 
single vector field. 

Proposition 3 

For a system of arbitrary number of interacting scalar fields minimally coupled to gravity 
and described by nontrivial Lagrange equations of motion, the full stress tensor cannot 
vanish for all solutions if it is not identically zero. 

Proof. 

For minimally coupled scalars any Lagrangian is of the form L(<f) a , <p a ^-, g^v) and is free 
of the metric connection and this results in the stress tensor free of the second order 
covariant derivatives, i.e. T^ u = T^ u ((p a , QiAt , g^ u ). This means that on a Cauchy surface 
the stress is determined solely by the initial data and thus can be given (since it 
does not vanish identically) any prescribed value (apart from some restrictions) inde- 
pendently of the field equations. By continuity, the stress tensor is also different from 
zero in some neighbourhood of the Cauchy surface. 

In the case of nonminimal coupling the proposition does not apply directly. However, all 
known cases comprise the scalar field of scalar-tensor gravity theories (all generalizations 
of Brans-Dicke theory), the conformally invariant scalar field (which mathematically is 
merely a special case of scalar-tensor gravity) and a scalar field arising in restricted 
metric nonlinear gravity theories (Lagrangian being a smooth function of the curvature 
scalar) via a suitable Legendre transformation. The nonminimal coupling takes the form 
f(4>)R, the Lagrangian may or may not contain a kinetic term for 0. Then one intro- 
duces a new spacetime metric by means of a Legendre map (actually this map is a 
conformal map of the original metric g^) and suitably redefines the scalar, cf) = <f)(4>). 
The mapping is commonly denoted as a transition from Jordan frame, i.e. the system 
{gfj,u,<P), to Einstein frame consisting of g^ v and cf). In Einstein frame <fi is minimally 
coupled to g^y (which should be regarded as the physical spacetime metric) [HI E] and 
Proposition 3 works. 

Proposition 4 

For a single vector field in Minkowski spacetime there exist solutions in the space of 
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all solutions vanishing sufficiently quickly at spatial infinity for which T^ V {A) ^ 0. 



Proof. 



A first order Lagrangian L^A^, A^) in Cartesian coordinates generates T lJLV {A a ^ A a ^, A a ^) 
explicitly depending on second derivatives except the case where L depends on A^ v only 
via F^y. Let S be a Cauchy surface with a unit timelike normal vector n u , n u n u = — 1. 
The initial data for A^ form a set Ca = {A^, n u A^ u } of functions given on S. On the 
boundary 2-sphere dS at spatial infinity the initial data vanish. 

One integrates the Belinfante-Rosenfeld identity (eq. (A. 12) in Appendix), expressed 
in Cartesian coordinates, over S, 



On S the tensor K^ ua being a linear combination of the classical spin tensor S^ va , eqs. 
(A.8-A.9), is determined by the initial data, K^ va = K^ voi {Ca)- Then applying the 
antisymmetry K^ VOL = —K mv and the Gauss' formula one can replace the integral of 
n u d a K^ ua over S by the integral of K^ va over dS and the latter is zero. Let A^ be a 
solution of E u = corresponding to the initial data Ca- Eq. (47) reduces to 



The canonical energy-momentum tensor fr* u , eq. (A. 7), is also determined on S by the 
initial data, whereas the stress tensor is determined by the solution A^. Since t^ v is not 
zero in general, one can choose such initial data Ca that the integral of ^"(Ca) does not 
vanish. Then the conserved total 4-momentum of the field for the solution is different 
from zero, 



The Proposition is proved. 

In this form the proof works only in Minkowski spacetime. Let A 0fM be a solution in 
flat spacetime for which eq. (49) holds. Making a small perturbation A^ = A 0fl + e M and 
9ni> = Vnu + hfu, (where may be a solution to Einstein field equations with some mat- 
ter source, not necessarily equal to A^) one finds that is altered by a small quantity 
T^ v which is linear in e M and h^. One then concludes that for spacetimes which are 
sufficiently close to flat space and for solutions being a small perturbation of A ^, the 
stress tensor is nonzero too. 

Furthermore, the proof of the Proposition may be generalized to at least one class of 
curved spacetimes. Let a spacetime metric admit a covariantly constant vector field 




(47) 




(48) 




(49) 
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kfM-u = 0. For example, if the spacetime is empty, R^ u = (and d — 4), then k^ is null 
and the spacetime represents the plane-fronted gravitational wave. In such a spacetime 
one considers a Cauchy surface extending to the spatial infinity. As previously, the initial 
data for are Ca = {A^, nUj ^^} with n u the unit normal vector to S; the data vanish 
on the boundary dS at spatial infinity. The tensors K^ va and are determined on S 
by the initial data whereas T^ v can be evaluated only for the solution corresponding to 
C A . 

One multiplies the identity (A. 12) by k^ and integrates it over 5* for a solution. The 
integral of the divergence term, k^nJSI a K lLVa = n v V a (k^K^ ua ^) , is equal to the integral 
over dS and hence is zero. One then arrives at an expression analogous to (48), 

[ k^T^ -t" u )n u dS = 0. (50) 
J s 

Once again, choosing such initial data that the integral of k^ v n v is nonzero, one gets 
T tlv {A) ^ 0. Unfortunately, the class of spacetimes admitting a covariantly constant 
Killing field is rather narrow. 

Finally we remark that Propositions 3 and 4 cannot be further strengthened, i.e. one 
cannot exclude the situation that the stress tensor vanishes in the whole spacetime for 
some particular solutions of the Lagrange equations of motion. In fact, at least one 
counterexample is known: a nonlinear massive spin-two field generated by a higher 
derivative gravity theory. For this field T^ u = in the spacetime of a plane-fronted 
gravitational wave [TU] . 

4 The canonical energy— momentum tensor in a 
curved spacetime 

In many theoretical investigations in classical and quantum field theory in Minkowski 
spacetime one employs the canonical energy-momentum tensor due to its conceptual 
simplicity and "naturalness". Yet the variational stress tensor is constructed in a way 
which is quite artificial in flat spacetime and one must appeal to arguments from outside 
the Lagrange formalism to show that the latter rather than the former is the physical 
energy-momentum tensor |H . 

The canonical energy-momentum tensor for a field ipA is 

t;^)E^-^— . (si) 

Taking divergence of the tensor and employing Ricci identity (the formula after eq. 
(A. 15)) and eq. (A. 14) one finds that 

dL 

V v yM = + 1 r—Z/ a R a ^. (52) 
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This means that even if the field equations hold, E A = 0, the canonical tensor is not 
conserved in a curved spacetime. This is the case of the electromagnetic (Maxwell) field, 



V„y(A) = - A a F R a 



(53) 



The above and other bizarre properties of the canonical tensor are due to the fact that 
the tensor, as defined in eq. (51), does not fit well to the variational formalism of field 
theory. In fact, any full divergence term in the field Lagrangian does not contribute to 
E A nor to T^iip), thus whether one discards such terms or not whilst evaluating these 
quantities does not affect the final outcome. Yet divergence terms in L do contribute 
to the canonical tensor. For example, let for a vector field L = V a ( i /'(A / M A[ ) J 4 Q ), then 



This example illustrates a generic feature of t^ u for vector fields: if the Lagrangian is a 
total divergence, then T 11V {A) = = and the Belinfante-Rosenfeld identity (A. 12) 
yields 



with a nonvanishing spin tensor. This vector field has no physical propagation and 
carries no energy, yet t^ u apparently attributes to it some nontrivial feature. This 
fact convincingly shows that the canonical energy-momentum tensor is not a physical 
quantity. 

Appendix 

For the reader's convenience we rederive here the Belinfante-Rosenfeld identity for any 
vector field generalized to a curved spacetime since it cannot be easily found in the 
literature in the form suitable for the present paper. The identity relates the varia- 
tional and the canonical energy-momentum tensors with the classical spin tensor and 
the Lagrange equations. Our approach is closest to that in [HI E21 HB] , f° r a different 
approach see [TU HH1 HE! • (The Belinfante-Rosenfeld identity for the vacuum electro- 
magnetic field, i.e. Maxwell's equations hold, in a curved spacetime is implicitly given 
in |17j.) We also recall the derivation of the Noether identity for any matter tensor 
field arising from the diffeomorphism invariance of the field action in the case of any 
metric theory of gravity. For an arbitrary vector matter field we show the explicit rela- 
tionship between the generalized Belinfante-Rosenfeld identity and the Noether identity. 

Let a classical vector matter field on a spacetime (M,g^ u ) be described by a 
scalar Lagrangian L. The Lagrangian may be written either in the explicitly covariant 



t^(A) = V a (fg» u A a - fg m A v ). 




(54) 
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way as a function of tensors or in the noncovariant form as a function of tensors and 
their partial derivatives, 

L(A„ A^, <f#) = L(A„ - n^Ar, g^) = L\A„ A^, g a \ g a \ v ). (A.l) 

The Lie derivative of the scalar density (of the weight +1) \J—gL with respect to an 
arbitrary vector field £ M may be evaluated either with the aid of the formula 



-gL) 



dA a 



dA a , 







(A.2) 



or by Lie differentiating it as a composite function, 



dA,. v 



OA. 



dg^, v 



dA 



dg afS 



-d v 



dA^ 



Oil 



(A.3) 



One uses the noncovariant form L' in order to arrive at the variational stress tensor and 
in the other differentiations one replaces L' by L. By applying C^A^ = ^ A^ v + A u ^ u ;ll 
and C^g a/3 = — 2£ ( - a;/3 ) one recasts (A.3) to a form explicitly exhibiting the dependence 
of the Lie derivative on £ M and its first and second order derivatives. First, denoting 

= and E*[L(A)] = £--vJ-^ 



dA a .p 

the Euler operator, one finds 

d(y/=gL') / OL 



dA„ 



OA 



dA,, 



+ 



dL OA, 



a;0 



and 



pAp dA a .p dA^ 
(d{y/=gL')\ 



, f dL 



dA~ d r <* 



dA„ v 



-g{d^, v -d^Y%\ 



Then the first square bracket in (A.3) equals ^J—gE^ 1 . The second square bracket is just 
— \\/~ 9 T a p{A). The first term in the third square bracket is also simple, 

d(y/=gL') _ ^— dL 



dA^ 



OA,./ 
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whereas the second term in this bracket is a complicated tensor density. After several 
manipulations with the last bracket and upon equating (A. 2) to (A. 3) one gets a scalar 
identity of the form 

+ ZkvQ 1 " + Up»)R a{ ^ = 0- (A.4) 

The identity holds for any vector fields and At any spacetime point the tensors 
£ M; „ and are independent quantities and therefore their coefficients must identically 

vanish. Thus the identity splits into a cascade of 4 + 16 + 40 identities. The identities 
Ra(v) = o are purely algebraic functions of of A^, g^ u and d^ v and are completely trivial 
in the sense that they hold for any values of these 3 tensors. The identities P M = are 
also algebraic functions of these tensors and R a p^ v (arising from Ca;^]) an d are trivially 
satisfied provided d^ u equals dL/dA^. u . Only Q^ u = are nontrivial and read 

where 



and = and N^" = d^ v \ One introduces the Pauli's canonical energy- 

momentum tensor 

tS(A) = d^L - A a ,,— — (A.7) 
and splits Kv va into its symmetric and antisymmetric parts w.r.t. to the first two indices, 

tf[H° = \A»d va - A v dT) = -S^ ua , (A.8) 



j^-(fj,v)a _ A a ^ v ) _ d a( - M A u ^ = S a( - pv \ (A. 9) 

The tensor £V^q may be interpreted as the density of the classical spin (helicity) of the 
vector field. In fact, for the electromagnetic field in vacuum in Minkowski spacetime the 
spin density tensor is defined as 

For L-^yi = — l/{lQn)F a pF a P and assuming that Maxwell equations hold, F alS $ = 
0, and the field vanishes sufficiently quickly at spatial infinity, one can remove a full 
divergence term from the spin density tensor since it does not contribute to the total 
spin of the electromagnetic field. As a result, 

S!Z° = —(A^F ua - A v F m ) = A»d va - A u d» a , (A. 11) 
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[IHlini. The tensor K^ a is antisymmetric, K^ a = -K^ au . Finally (A.5) takes on the 
form 

QT = - tr + A»E V + ^V a (S^ a + S a ^ + S av ") = 0. (A. 12) 

These are the Belinfante-Rosenfeld identities fHl EDj for any vector field, generalized 
to a curved spacetime. Notice that the identities arise due to the fact that the field 
Lagrangian is a scalar, hence the action is invariant under infinitesimal coordinate trans- 
formations (with appropriate boundary conditions). 

The proof of Propositions 1 and 2 in sect. 2 is based on the famous second Noether 
theorem (also named "the Noether identities" , "the strong conservation laws" or "Bianchi 
identities for matter") [2U 122 EHj ; our approach is based on [23 EH- Let ipA denote an 
arbitrary tensor matter field (or a set of tensor fields) with a collective index A, described 
by a scalar Lagrangian L(ip) = L(ip a, ip a-,^ g^v, R a /3^u), i-e. one admits a possible non- 
minimal coupling to the curvature. The action integral for for Einstein gravity (actually 
one may envisage any metric theory of gravity since the derivation goes then without 
alteration) and the matter field is (we set c = 87rG = 1) 

S= [ (-R + L)^d 4 x. (A.13) 
Jn 2 

The Euler operator for the equations of motion for ipA is 

The action (A.13) is diffeomorphism invariant and in particular it remains unchanged by 
an infinitesimal point transformation x ,fl = x M + ^(x) with arbitrary £ M which vanishes 
on the boundary of the domain Q. The integration domain is thus mapped onto itself 
whereas the variations of the fields are given by Lie derivatives, 5ipA = —C^a an d 
SgW — —£^gP v These variations can be expressed in terms of coefficients Z A aiip), 
being linear functions of i[) A , 

CtfA = rVU;a + VaO, (A.15) 

which also appear in the covariant derivative, 
and in Ricci identity, 
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The invariance of the action implies 

= 5S= [ [\{G, V - T^SgT + E A 5^ A ]^g~d A x (A.16) 

plus a surface integral which vanishes due to the boundary conditions. Applying (A. 15), 
dropping again a total divergence and making use of the ordinary Bianchi identity 
G v = one finds that the integrand is of the form £ M B ^(ifj) • Vanishing of the integral 
and arbitrariness of £ M in the interior of Q entail _B M = or 

E A ^A-, a = V (Tj + E A Z A P a ). (A.17) 

This is the Noether identity for any classical tensor matter field, valid for any dimension 
d > 3. Notice that the second Noether theorem derived in refs. does not include 

the variational (stress) energy-momentum tensor T M1 , (?/>). The tensor is absent there 
because it is (implicitly) assumed in these papers that the matter action is invariant 
under a symmetry transformation of the matter field (a gauge transformation) alone 
and the metric is left unaltered. 

In the case of a vector field the coefficients are Zj 3 a = S^A a and the four Noether 
identities 

E»F atl = T a ^ + E%A a , (A.18) 

where F a p = A^ a — A a;/ 3, can be derived from the Belinfante-Rosenfeld identities (A. 12). 
In fact, a direct and quite uphill calculation proves that 

Q a ^ = V ~ F ^ Efl + AoE^ = 0. (A.19) 
Analogous relationships exist for tensor fields. 
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